arXiv:1507.00979vl [math.PR] 3 Jul 2015 


Bounds of the accuracy of the normal approximation to the 
distributions of random sums under relaxed moment conditions* 

V. Yu. Korolev! A. V. Dorofeeval 


Abstract: Bounds of the accuracy of the normal approximation to the distribution of a sum of 
independent random variables are improved under relaxed moment conditions, in particular, under the 
absence of moments of orders higher than the second. These results are extended to Poisson-binomial, 
binomial and Poisson random sums. Under the same conditions, bounds are obtained for the accuracy 
of the approximation of the distributions of mixed Poisson random sums by the corresponding limit law. 
In particular, these bounds are constructed for the accuracy of approximation of the distributions of 
geometric, negative binomial and Poisson-inverse gamma (Sichel) random sums by the Laplace, variance 
gamma and Student distributions, respectively. All absolute constants are written out explicitly. 
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1 Introduction 


1.1 The history of the problem and aims of the paper 


Let Xi,X 2 , ... be independent random variables with EAj = 0 and 0 < = <7^ < oo, i = 1,2,... For 

n E N denote Sn = Xi X^, = a\ + ... + Let <I>(x) be the standard normal distribution 

function. 


X 



— QO 


Denote 

Xn = supIP!^^ < xBn) - ^>(a:)|. 

X 

Let Q be the class of real functions g{x) of x E M such that 

• g{x) is even; 

• g{x) is nonnegative for all x and g{x) > 0 for x > 0; 

• g{x) and x/g{x) do not decrease for x > 0. 

In 1963 M. Katz |25j proved that, whatever 5 E ^ is, if the random variables Ai, A2 ,... are identically 
distributed and EXfg{Xi) < 00, then there exists a finite positive constant Ci such that 


An ^ Cl ■ 


EAfg(Ai) 

(Xi9{crii/n)' 


( 1 ) 
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In 1965 this result was generalized by V. V. Petrov [lOj to the case of non-identically distributed summands 
(also see | 11 | 1 : whatever <7 € ^ is, if ^XfgfyXi) < 00 , i = 1 ,..., n, then there exists a finite positive constant 
C 2 such that 




C 2 






2=1 


(2) 


Everywhere in what follows the symbol ]I(^) will denote the indicator function of an event A. For 
£ G ( 0 , 00 ) denote 


- /t. 'X 

^ ^Bn), Mnie) = < eBn). 

i=\ i=l 

In 1966 L. V. Osipov [9] proved that there exists a finite positive absolute constant C 3 such that for any 
£ G (0, 00 ) 

An ^ C 3 [Ln(£) + Mn{£)\ (3) 

(also see |12| . Chapt V, Sect. 3, theorem 7). This inequality is of special importance. Indeed, it is easy to 
see that 

n 

Mn(e) < eBn) ^ e. 

i=i 

Hence, from (3) it follows that for any e G (0, oo) 

An ^ C *3 (e + Ln{£)) ■ (4) 


But, as is well known, the Lindeberg condition 


lim Ln{£) = 0 for any e G (0,oo) 

n—>oo 


is a criterion of convergence in the central limit theorem. Therefore, in terminology proposed by 
V. M. Zolotarev [35], bound (4) is natural, since it relates the convergence criterion with the convergence 
rate and its heft-hand and right-hand sides converge to zero or diverge simultaneously. 

In 1968 inequality (3) in a somewhat more general form was re-proved by W. Feller |22j . who used the 
method of characteristic functions to show that C 3 ^ 6 . 

A special case of (3) is the inequality 


A„^C^L„(1) + M„(1)]. 


In the book m it was demonstrated that C 3 ^ ‘^C^. 

For identically distributed summands inequality (5) takes the form 




An ^ min < 1, 

O'a 


1^11 


Oly/n 


(5) 


( 6 ) 


In the papers [291130] L. Paditz showed that the constant can be bounded as < 4.77. In 1986 
in the paper m he noted that with the account of lemma 12.2 from | 2 |, using the technique developed 
in |291130| . the upper bound for C 4 can be lowered to 6*4 < 3.51. 

In 1984 A. Barbour and P. Hall |18j proved inequality (5) by Stein’s method and, citing Feller’s result 
mentioned above, stated that the method they used gave only the bound C'^ ^ 18 (although the paper 
itself contains only the proof of the bound C 3 ^ 22). In 2001 L. Chen and K. Shao published the paper 
m containing no references to Paditz’ papers ESI EDI El] in which the proved inequality (5) by Stein’s 
method with the absolute constant C'^ = 4.1. 

In 2011 V. Yu. Korolev and S.V. Popov |26j showed that there exist universal constants Ci and C 2 
which do not depend on a particular form oi g ^ Q, such that inequalities (1), (2), (5) and ( 6 ) are valid 
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with Cl = (74 ^ 3.0466 and C 2 = C 3 ^ 3.1905. This result was later improved by the same authors in 
the papers PE], where it was shown that Ci = C 2 = C 4 = C 3 ^ 2.011. 

Moreover, in the paper |7] lower bounds were established for the universal constants Ci and C 2 . 
Namely, let g be an arbitrary function from the class Q. Denote by Tig the set of all random variables X 
satisfying the condition EX'^g{X) < 00 . Denote 


C* =sup sup 
g&G 

1=1,... ,n 


XnBlgjBn) 

ELi ^X‘fg{Xiy 


It is easily seen that C* is the least possible value of the absolute constant C 2 that provides the validity 
of inequality (2) for all functions g ^ Q at once. In the paper [7| it was proved that 


C* 


^ sup 
^>0 


1 


1 + z^ 




= 0.54093. 


The aim of the present paper is to improve and extend the results mentioned above. First, we will show 
that one can take C 3 = C 3 . Second, we will sharpen the upper bounds of the absolute constants mentioned 
above. Third, we will extend these results to Poisson-binomial, binomial and Poisson random sums. Under 
the same conditions, bounds will be obtained for the accuracy of the approximation of the distributions of 
mixed Poisson random sums by the corresponding limit law. In particular, we will construct these bounds 
for the accuracy of approximation of the distributions of geometric, negative binomial and Poisson-inverse 
gamma (Sichel) random sums by the Laplace, variance gamma and Student distributions, respectively. 
All absolute constants will be written out explicitly. 

Along with purely theoretical motivation to sharpen and generalize known results, there is a somewhat 
practical interest in the problems considered below. Poisson-binomial, binomial and mixed Poisson (hrst 
of all, geometric) random sums are widely used as stopped-random-walk models in many fields such 
as financial mathematics (Cox-Ross-Rubinstein binomial random walk model for option pricing [20)1. 
insurance (Poisson random sums as total claim size in dynamic collective risk models |21| . binomial random 
sums as total claim size in static portfolio risk models, geometric sums in the Pollaczek-Khinchin-Beekman 
representation of the ruin probability within the framework of the classical risk process Ell), reliability 
theory for modeling rare events |24) . It is now a tradition to admit that the distributions of elementary 
jumps of these random walks may have very heavy tails. The problems considered in the present paper 
correspond to the situation where the tails may be as heavy as possible for the normal approximation 
to be still adequate. Moreover, the bounds obtained in this paper partly give an answer to the questions 
how heavy these tails can be for the normal approximation (or scale-mixed normal approximation) to be 
reasonable. 

The paper is organized as follows. In Section 1.1 we prove that in inequalities (l)-(5) the absolute 
constants coincide and that the values of these constants are determined by that of In Section 2 the 
upper bound of is sharpened. In Section 3 the analogs of inequalities (1), (2), (3) and ( 6 ) are proved 
for Poisson-binomial and binomial random sums. In Section 4 the results obtained in Section 3 are used to 
construct the analogs of (1) and ( 6 ) for Poisson random sums. The results of Section 4 are used in Section 
5 to obtain bounds for the accuracy of the approximation of the distributions of mixed Poisson random 
sums by the corresponding limit law. In particular, here these bounds are constructed for the accuracy of 
approximation of the distributions of geometric, negative binomial and Poisson-inverse gamma random 
sums by the Laplace, variance gamma and Student distributions, respectively. 


1.2 On the coincidence of the absolute constants in inequalities (1)—(5) 

The main result of this section is the following statement. 

Lemma 1. For any e € (0, 00 ) 


-bn(l) + Af„(l) ^ Ln{s) -|- Mn{s). 


(7) 
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Proof. For e = 1 the statement is trivial. Let e < 1. Then 


Ln{^) + -Lfn(l) — L„(e) + M„(e)+ 

+ ^ Y, \Xj?l{eBn < |X,| <Bn)-^Y ^ l^il < ^n)- 

"■ i=i " i=i 

But 

^ n 1 ^ 

^ \Xj\ < Bn) — Y^ ^XjI{eBn ^ \Xj\ < Bn) ^ 

n j=i n 

1 n 1 ^ 

^ -4 ^ l^il <^n)-J^Y ^ |X, | < = 0, 

n j=i n 

therefore, in the case e < 1 inequality (7) is proved. 

Let now £ > 1. Then 

-Ln(l) + Mn{l) = Lnis) + Mn{£) + 

+ ^ Y ^^PiBn ^ l^il < sBn) -^Y l^jPiBn ^ |X,| < sBn). 

"■ i=i "■ i=i 

But 

^ n 1 ^ 

Y ^Xfl{Bn ^ |Xj| < eB„) --^Y \XjP{Bn ^ \Xj\ < eBn) ^ 

n j=i n 

1 n 1 ^ 

^ -4 ^ l^il < ^Bn) -JhY ^Y^(Bn ^ |X, | < eBn) = 0, 

n j=i n 

that is, the statement of the lemma holds for e > 1 . as well. 

Corollary 1. The absolute constants in inequalities (3), (4), (5) and (6) can be taken identical, that 
is, if inequality (5) holds with C'^ ^ Cq, then inequalities (3), (4) and (6) hold with C 3 ^ Cq n (74 ^ Cq. 

Remark 1. In the paper [7] it was shown that if inequality (5) holds with (7^ ^ (7o, then inequalities 
(1) and ( 2 ) hold with Ci ^ Cq, i = 1,2. 

So, in the evaluation of the constants in the above inequalities, the constant C'^ in inequality (5) plays 
the determining role: if a particular upper bound C 3 ^ Cq is known, then in all the rest inequalities 
(l)-(4) and ( 6 ) one can let Ci ^ Cq, i = 1,2, 3,4. That is the reason for us to focus on sharpening the 
upper bound for (7^. 


2 Sharpening of the upper bound for the constant C'^ 

2.1 Auxiliary results 

For X ^ 0, n G N and i = 1,..., n denote 

n 

Y,{x) = Bn^Xil{\X,\<{l + x)Bn), Y, = Y,iO), Wn(x) = Y^ii^)^ H"n = lLn(0). 

i=l 


Since EAj = 0, we have 


EAil(|A,| < (l + x)R„)| = |EA,]I(|A,| ^ (l + x)R„)|. 


( 8 ) 
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By the definition of the random variables Yi{x) the relation 


i=l 


1 

^ i=l 


= 1 


(9) 


holds. Denote 


K = 


17 + 7v^ 

27 


< 1.3156. 


Lemma 2. 1°. For any n G N, x ^ 0 and p G [1, K] there holds the inequality 

^ E|yi(x) - Ey,(x)|3 < min {iLM„(l + x), pM„(l + x) + (5 " P) ^n(l + x) 
i=i ^ 

2°. For any n G N and x ^ 0 t/zere hold the inequalities 

1 - 2F„(1 + x) ^ DlL'n(x) ^ 1. 

3°. Let Mn{l) = 7L„(1), 7^0. Then for any n G N there holds the inequality 


Y, E|Ei - Ey,|3 < Ln{l) min { 7 ^ 7 , 7 + 4} . 

i=l 

The PROOF based on the results of [231 M and [ 8 ] was given in | 6 ] . 
Lemma 3. 1°. Let q > 0. Then 


sup | 4 >(( 7 x) — 4>(x)| = 


$ q 


I Ingr^ 

— 1 


- 4> 


' lny2 

— 1 






2°. Let a G M. Then 


fq-l)lnq 
7 r{q + 1 ) 


exp 


{ - mind, max {,, L - 1) 


\/27re 


1 - 

q- 


sup |<h(x + o) — <h(x)| = 2 <hf—'j — 1 ^ 
X V 2 / 


v^’ 


The elementary PROOF of this lemma is based on the Lagrange formula and the easily verihable fact: 
if F(x) and G(x) are two differentiable distribution functions, then sup |F(x) — G'(x)| is attained at those 

X 

points X, where F'(x) = G'{x) (also see [H], p. 143). 

Lemma 4. Assume that L„(l) ^ A for some A G (0, ^). Let 

2 


B{A) = 


(1 + VI - 2.4)VI - 2 A' 


Then 


1 ^ 


1 


V 1 + B{A)Ln{l). 


For the PROOF see [ 6 ]. 

Lemma 5. Let X be a random variable with EX^ < 00 . Then 


sup 


X-EX . , 


V sup 
^>0 


1 


l + z^ 


-4>(-V 


= 0.54093. 


For the PROOF see, e. g., the book |2] and the papers mm- 
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2.2 General case 


Theorem 1. Let n € N, the random variables Xi,... be independent, EXj = 0 and 0 < EX? < oo, 
i = 1,... ,n. Let 7 = M„(l)/L„(l). Then there exists a finite positive number 6 * 1 ( 7 ) depending only on 7 
such that 

An ^ (1 + 7 ) 6*1 ( 7 )Ln(l). 

Moreover, the upper bounds for 61 ( 7 ) are presented in table 1. 


7 

61 ( 7 ) + 

7 

61 ( 7 ) + 

7 

61 ( 7 ) ^ 

7^0 

1.8627 

7 + 1 

1.5605 

7+10 

0.9393 

7 ^ 0.1 

1.8587 

7 + 2 

1.3488 

7 ^ 100 

0.6067 

7 ^ 0.5 

1.7244 

7 + 5 

1.0836 

7 —>■ 00 

0.5583 


Table 1: Upper bounds for 61 ( 7 ). 

Corollary 2. Under the eonditions of theorem 1, inequalities (2) — (5) hold with 62 = 63 = 6*3 ^ 
1.8627. 

Proof of theorem 1. For any y € M the event {S*n < yBn} implies the event 


{Wn <y}U {|Xi| Bn} U...U {\Xn\ ^ Bn}, 
whereas the event {Wn < y} implies the event 


{5n < yBn} U {|Xi| ^Bn}U...U {|Xn| ^ Bn}. 


Therefore, 

Hence, for any y G 
where 


sup |P(Sn < yBn) - P{Wn < y)K P(|Xi| ^ Bn). 


2 = 1 


Qi = sup 
y 

Q 2 = sup 
y 


An ^ Ql + Q2 + Q3, 


P J Wn- EWn ^ y-PWn '^ _ ^ f y - 


( 10 ) 


$ 


VWv, 

y-EW, 

\/DW^ 






-^{y) 


Q3 = ^P{\Xi\fi Bn). 


2=1 


Consider Qi. By virtue of the Berry-Esseen inequality with the best known upper bound of the absolute 
constant (see |15) 1 we have 

0.5583 


Ql ^ 


{DWnf‘' .=1 


^ElXi-ETil^. 


Assume that Lri(l) ^ < i. Then in accordance with statements 2° and 3° of lemma 2 


Ql ^ 


0.5583 • min{X 7 , 7 + 4} Ln(l) 

(1 - 2A)3/2 ■ 


( 11 ) 


Consider Q 2 . We obviously have 


Q 2 = sup 
y 


y-PWn \ 

^/DWn ) 


$ (y - EWn) + ^ ( 2 / 


EWn) - 4>(y) 
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^ sup 
y 

= sup <I> ( 
y ' 


y-^Wn 

\/DT^ 


- $ (y - EWn) + sup I - EWn) - ^>(y)| = 
y 


- $ (y) + sup I ^>(y - EWn) - ^(2/)| = Q 21 + Q22- 


According to statement 2° of lemma 2, DWn ^ 1. Therefore, by virtue of statement 1° of lemma 3 and 
lemma 4, there holds the inequality 


Q 21 ^ 




-1 < 


_ 2 T^( 1 ) _ 

/27re(l - 2A)(1 + y/1 - 2A)' 


Consider Q 22 - By virtue of (8) we have 


n . n 1 ^ 

\^Wn\ = < Bn)\ = |EX,I(|Ai| ^ Bn)\ ^ 


- ft ^ ft 

^ Bn) = Lnil). 


Therefore, by statements 2° of lemma 2 and 2° of lemma 3, 


Q 22 ^ 


Tn(l) 


From (12) and (13) it follows that 


Tn(l) 


/e(l - 2 A )(1 + Vl - 2 A) 


Finally, by the Markov inequality 


Q3 = Y1 P(I^*I ^Bn)^-^Y1 ^ Bn) = Lnil). 

i=l i=l 


So, from (10), (11), (14) and (15) we obtain 


^ Tjj(l) 1 + 


/e(l - 2 A )(1 + VI - 2 A) 


Introduce the function 


H,ij,A) = l + - 


V(1 - 2A)(1 + VI - 2A) 


0.5583 • minjiFy, 7 + 4} 
(1 - 2A)3/2 


0.5583 • minjiFy, 7 + 4} 
(1 - 2A)3/2 ■ 


For any 0 V A < ^ we have the inequality 

An ^ L„(l) • max ji^Vq, A), ^ 

This follows from (16) if Lnil) ^ A and from lemma 5 otherwise. 
Now, with the account of the equality 


Lnil) = 


iLnil) + 'yLnil)) _ Lnil) + Af„(l) 
1+7 1+7 
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we have 


C*i(7) ^ min max 
o^A<i 


0.541 ) 

1 + 7 ’ A{1 + 7 ) / ’ 


The computation by this formula yield the values presented in table 1. Note that the hrst function of A 
inside the minimax is increasing whereas the second one is decreasing. Hence, the value of the minimax 
is delivered by the unique solution of the equation 


Hi{'j,A) 0.541 
1 + 7 H(1+7)' 


For 7 > 13 (we have 7 + 4 < K'j) both functions decrease in 7 , that is, the minimax value decreases. 
Therefore, the corresponding part of table 1 is obtained by the evaluation of the bound for Ci{'y) at one 
point. The part of table 1 corresponding to 0 ^ 7 ^ 13 is obtained by numerical optimization of a finite 
interval. The theorem is proved. 


2.3 Special cases 

Using the best current upper bound Cq ^ 0.4690 for the absolute constant in the Berry-Esseen inequality 
for identically distributed summands (see m), the following statement can be obtained in the way similar 
to the proof of theorem 1 . 

Theorem 2. In addition to the assumptions of theorem 1, let the random variables Xi,X 2 ,... be 
identically distributed. Then there exists a finite positive number ( 72 ( 7 ) depending only on 7 such that 

An ^ (1 + 7)C'2(7)Tn(l)- 

Moreover, the upper bounds for C 2 {'y) are presented in table 2. 


7 

62 ( 7 ) < 

7 

62 ( 7 ) ^ 

7 

62 ( 7 ) < 

7^0 

1.8546 

7 + 1 

1.4793 

7+10 

0.8292 

7 ^ 0.1 

1.8338 

7 + 2 

1.2540 

7 + 100 

0.5147 

7 > 0.5 

1.6608 

7 + 5 

0.9781 

7 —>■ c» 

0.4690 


Table 2: Upper bounds for 6 * 2 ( 7 ). 


Proof. Using the reasoning similar to that used to prove theorem 1, it is easy to see that 


where 


H2{'y,A) = 1 + 


. fH2il,A) 0.541 ) 

07A<i 11+7 ^(1 + 7)/ 

2 \ 0.4690 • minjiFy, 7 + 4 } 




1 + - 


(1 - 2H)3/2 


^e(l- 2 H)(l + VI - 2 A) J 
The computations by these formula yield the values of the upper bounds for 62 ( 7 ) presented in table 2. 
The theorem is proved. 

Corollary 3. Under conditions of theorem 2, inequalities (1) and ( 6 ) hold with 61 = 64 V 1.8546. 
Theorem 3. In addition to the conditions of theorem 1, let the random variables Xi,X 2 , ■ ■ ■ have 
symmetric distributions. Then there exists a finite positive number 63 ( 7 ) depending only on 7 such that 


An ^ (1 + 7)63(7)Ln(l)- 


Moreover, the upper bounds for 63 ( 7 ) are presented in table 3. 
























7 

63 ( 7 ) < 

7 

63 ( 7 ) < 

7 

63 ( 7 ) ^ 

7 + 0 

1.5769 

7 + 1 

1.3033 

7+10 

0.7433 

7 + 0.1 

1.5749 

7 + 2 

1.1115 

7 + 100 

0.5808 

7 + 0.5 

1.4532 

7 + 5 

0.8729 

7 —>■ 00 

0.5583 


Table 3: Upper bounds of 6*3 ( 7 ). 


Corollary 4. Under the conditions of theorem 3, inequalities (2) — (5) hold with C 2 = C 3 = C'^ ^ 
1.5769. 

The PROOF of theorem 3. In the case under consideration instead ( 8 ) we have 

^ ^ 0.5583M„(1) 

" (1 - 2A)3/2 ’ 

and <522 = 0, since EWn = 0. Therefore, the bound 

2 


^ Ln{l) ( 1 + 


+ 


holds. Thus, 


where 


0.5583Mn(l) 
V^27re(l - 2^)(1 + Vl - 2^)7 ' (1 - 2A)3/2 

0.541 


^* 3 ( 7 ) ^ min max 
07A<i 


1 + 7 ’ yl(l + 7 ) 


= 1 + 


+ 


0.55837 


V^27re(l - 2A){1 + Vl - 2A) (1 - 2yl)3/2 ' 

The computations by the above formulas yield the values of the upper bounds for 6 * 3 ( 7 ) presented in 
table 3. The theorem is proved. 

Theorem 4. In addition to the conditions of 3, let the random variables Xi, X 2 , ■ ■ ■ be identically 
distributed. Then there exists a finite positive number 64 ( 7 ) depending only on 7 such that 

An ^ (1 + 7)C'4(7)Ln(l)- 

Moreover, the upper bounds for 6 * 4 ( 7 ) presented in table 4. 


7 

64 ( 7 ) < 

7 

64 ( 7 ) < 

7 

64 ( 7 ) ^ 

7 + 0 

1.5645 

7 + 1 

1.2388 

7+10 

0.6591 

7 + 0.1 

1.5534 

7 + 2 

1.0373 

7 + 100 

0.4923 

7 + 0.5 

1.4018 

7 + 5 

0.7915 

7 —>■ 00 

0.4690 


Table 4: Upper bound for 64 ( 7 ). 


Corollary 5. Under the conditions of theorem 4, inequalities (1) and ( 6 ) hold with 61 = 6*4 ^ 1.5645. 
Proof of theorem 4. In the case under consideration 


64 ( 7 ) ^ min max 
07A<i 


774 ( 7 , 2 !) 0.541 

1 + 7 ’ yl(l + 7 ) 


where 


774 ( 7 ,^) 


1 + 


2 

V^27re(l - 22 l)(l + y/T^^) 


0.46907 
(1 - 2 ^) 3/2 ■ 


The computations by the above formulas yield the values of the upper bounds for 64 ( 7 ) presented in 
table 4. The theorem is proved. 
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3 The accuracy of the normal approximation to the distributions of 
Poisson-binomial random sums 

From this point on let Xi,X 2 , ... be independent identically distributed random variables with EXi = 0 and 
0 < EXf = < oo. Let pj G (0,1] be arbitrary numbers, j = 1,2, .... For n G N denote On = pi + . ■ ■+Pni 

Pn — (Pi) ■ ■ ■ jPn)- The distribution of the random variable 

^n,p^ = Cl T . . . + Cn; 

where Cl; • • •; Cn are independent random variables such that 

{ 1 with probability pj, 

, j = 1 ,... ,n, 

0 with probability 1 — pj, 

is usually called Poisson-binomial distribution with parameters n',Pn- Assume that for each n G N the 
random variables Nn,p^, Xi, X 2 , ■ ■ ■ are jointly independent. The main objects considered in this section 
are Poisson-binomial random sums of the form 


AAr„ . 

As this is so, if Nn » = 0, then we assume Sm =0. 

For J G N introduce the random variables Xj by setting 





with probability pj, 
with probability 1 — Pj- 


If the common distribution function of the random variables Xj is denoted F{x) and the distribution 
function with a single unit jump at zero is denoted Eq{x), then, as is easily seen, 

P{Xj < x) =pjF{x) (1 —pj)Eo{x), X G M, j G N. 


It is obvious that EXj = 0, 

E)Xj = ex] = pja^. 


In what follows the symbol = will denote coincidence of distributions. 
Lemma 6. For any n G N and pj G (0,1] 


(18) 


=Xi + ... + A„, (19) 

where the random variables on the right-hand side of (19) are independent. 

Proof. The characteristic functions of the left-hand and right-hand sides of (19) have the following 
forms 


= '^PXi+...+Xkit)PiNn,p^ = k) and (t) = Yl\pjPXj{t) + (I - Pj)]- 

k=0 j=l 

It suffices to make sure that the characteristic functions of the left-hand and right-hand sides of (19) 
coincide. 

We will use the method of mathematical induction. Basis: n = 1. 

PlPXiit) (1 -pi) =piipxi{t) (1 -pi). 
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Inductive step: we show that if the characteristic functions of the left-hand and right-hand sides of 
(19) coincide with n = m, then they also coincide with n = m -|- 1. 


m+l 


n + (1 - Pi)] = + {1 - Pj)]){pm+l^Xrr,+i{t) + (1 -Pm-rl)) = 

i=i i=i 

m m 

(1 ~ Pm+l) 'y ]] PXi-r...-rXfc(^)P(-^m,p„ = k) + Pm+lPXm+i (^) PXi-r...-rXfc (^)P(-^m,p^ = fc) 

/ i ;=0 / c =0 

m 

= (1 -Pm+i) ^V?Xi-r...-rXfc(t)P(iVm,p^ = /c) -hpi.. .pm+iv?Xi-r...-rx^+i + 


Ai=0 


m—1 


-Fp™+i(/?x^+i(t) ^ (^Xi-r...-rXfc(t)P(iVm,p^ =k) = 
k=0 

= (note that (pxi (^) = • • • = PXm+i (^) transform the last term) = 

m 

= {I - Pm+l)'^PXi+...+Xk{t)^{^ni,p^ = k) +pi . ..pm+WXi+.. 


k=0 


+Pm+l'^‘fXi+...+Xkit)P{Nm,p^ =k-l). 


k=0 


On the other hand, 


m+l 


‘fXi+...+Xkit)P{Nm+l,p^+l =k) = 

k=^ 

m 

'^PXi+...+Xk{t)P{Nm+l,p^^, = k) +Pi . ..pm+WXi+...+Xr.+i = 


Ai=0 

m 

~ ^ ^ V^Xi+...+Xfc (^)P({-^m,p^ = fc n ^m+1 = 0} U = k — 1 D <^m+l = 

k=0 

-\-pi . . .pm^iipXi-\-...-\-Xm+i = 
m m 

= [l - Pm+l)'^^Xi+...+Xk{t)P{Nm,p^ = k) + Pm+l^'PXi+...+Xk{t)P{N^,p^ = k - 1) + 
k=0 k=0 

+Pl ■ ■ ■Pm+lPXl+...+Xrr,+ l- 

Note that the right-hand sides of the above chain of equalities coincide. The lemma is proved. 
With the account of (18) and (19) it is easy to notice that 


= OnO- . 


( 20 ) 


Denote 


An,p„ = sup |P(5Ar„ ^^ < Xa^/K) - $(x) . 


Theorem 5. For any n G N and pj G (0,1], j G N, 


, 1.8627 ^,, 

^n,p„ ^ EX 
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Proof. By virtue of lemma 6 and relation (20) we have 


= sup |P(Xi + ... + < xa^/^) - d>(x)|, 

X 

and for the latter expression we can use the bound given in theorem 1 : 

sup |P(Xi + ... + Xn < xa\f 6 ^) - <I>(x)| ^ 

X 


^ 1.8627 


a‘^6^ 


^ n 1 ^ 


i=i 


= 1.8627 


a‘^9n 
= 1.8627 


^ n 1 ^ 


i=i 

a'^On 


” E|Xi|3i(|Xi| <u\/^) ” 


i=i 


i=i 


= 1.8627 


CJ 

1.8627 


^EXfldXil ^ a^n) + ^i=E|Xi|3l(|Xi| < a^n) 

^ V “n 


EXi minjcrv^, |^i|} = 


1.8627 


EXi m i n < 1 


^ 1^11 
’ (yy/K 


Q. E. D. 

Theorem 6. Under the conditions of theorem 5, whatever funetion g G G is such that ^X‘lg{Xi) < oo, 
there holds the inequality 

A„, < 1.8627 

" •y^iAasfTn) 

Proof. Let g be an arbitrary function from the class Q. With the account of the properties of a 
function gf G ^ it is easy to see that 

EXfl{\Xi\ ^ ay/K) = EX^^^IdXil ^ ay/K) ^ —^^EXfg{X,)I{\X,\ ^ ay/^) (21) 

9[^i) g{(yy/9n) 


and 


EXfl{\X^\ < ay/K) = EXfg{X^)^^Ii\X^\ < oyf^) ^ EXfg{Xi)I{\Xi\ < ayf^). (22) 

5 (^ 1 ) g{(yy/K) 

Substituting these estimates into the inequality 


An,p„ ^ 1.8627 


^EXfl{\X,\ ^ ay/^) + -^E\X^\H{\X^\ < ayfWn) 


obtained in the proof of theorem 5, we have 


(23) 


A 




< 


1.8627 

(y‘^g{(Tyfdn) 


[EXlg{Xi)l{\X^\ 




ay/Fn) + EXlg{Xi)l{\Xi\ < ay/Fn)] = 1.8627 


EXfg(Xi) 

(y‘^g{(yy/F) 


The theorem is proved. 

In particular, if pi = p 2 = • • • = P; then the Poisson-binomial distribution with parameters n G N and 
becomes the classical binomial distribution with parameters n and p: 




n,Pr 




P(iVn,p = A:) = C^p'^(l-p) 


n,p5 


\n—k 


/c = 0 , . . . 


n. 
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In this case 9n = np, so that ^Sn„p = npa‘^. Denote 

An,p = sup |P(S'Ar„_p < - $(x)|. 

X 


Estimates of the accuracy of the normal approximation to the distributions of binomial random sums 
(under traditional conditions of the existence of the third moments of summands) were considered in 
the paper |33| . where a conventional approach was used which is based on the direct application of the 
total probability formula and does not involve representation (19). Hence, in |33| estimates were obtained 
with the structure far from being optimal, containing unnecessary terms and unreasonably large values 
of absolute constants. 

Theorems 2 and 5 imply 

Corollary 6. For any n G N and p € (0,1] 






1.8546 




min < 1, 


l^il 

cr^/npj ■ 


Theorems 2 and 6 imply 

Corollary 7. Under the conditions of theorem 5, whatever function g £ G is such that ^X'lg{Xi) < 
oo, for any n G N and p G (0,1] there holds the inequality 


‘■n,Pr, 




1.8546 


a‘^g{a,/np) ’ 


4 The accuracy of the normal approximation to the distributions of 
Poisson random sums 

In addition to the notation introduced above, let A > 0 and Nx be the random variable with the Poisson 
distribution with parameter A: 

\ k 

p{Nx = k) = e-^ — , /sGNU{0}. 

k\ 

Assume that for each A > 0 the random variables Nx,Xi,X 2 , ■ ■ ■ are jointly independent. Consider the 
Poisson random sum 

Snx — Xi + ... + X]\[^. 

If Nx = 0, then we set = 0. It is easy to see that ESa = 0 and DS'a = Xa^. The accuracy of the 
normal approximation to the distributions of Poisson random sum was considered by many authors, see 
the historical surveys in [271 El. However, the authors are unaware of any analogs of the Katz-Osipov-type 
inequalities (1) and (6) under relaxed moment conditions. 

We will obtain a bound for 


Aa = sup IP [Sx < xaVx) — 4>(x)I. 

X 

For this purpose fix A and along with Nx consider the random variable Nn^p having the binomial 
distribution with arbitrary parameters n and p G (0,1] such that np = A. As this is so, the reasoning used 
above implies that 

DSnx = = fJ^A = a'^np. 

Therefore, by the triangle inequality, in accordance with corollary 6 we have 

Aa ^ An,p + sup |P(5Ar;, < x) - PiSN„,p < x)| ^ 


13 





1.8546 


^ x^^^^2min J 1, I+sup^Pf <xj|P(iV„,p = fc)-P(iVA = A:)| < 


a^/np 


1.8546 ^ . 

^ - 7 ^— EX^ mm 

cr^ 


k=o ^ j=i 

OO 


{ 1, ^ |P(A^n,p = k)- P(Xa = k)\. 


( 24 ) 


Estimate the second term on the right-hand side of (24) by the Prokhorov inequality |13) (also see |17) . 
p. 76), according to which 


|P(A^n,p = k) - P{Nx = k)\ ^ 2pmin{2, A}, 

k=0 

and obtain that for any n and p such that np = X, there holds the inequality 

Aa ^ min |l, + 2pmin{2, A}. 

Now, putting in (25) p = Xfn and letting n ^ oo, we obtain the final result: 
Theorem 7. For any A > 0 


Ax ^ 


1.8546 


(T^ 


EX^ min < 1 


Ml 

(tVX) 


(25) 


Using inequalities (21) - (23) to estimate A„^p in (24), we obtain the following result. 
Theorem 8 . Whatever function g £ Q is such that EX^g'(Xi) < oo, there holds the inequality 


Aa ^ 1.8546 


EXfg(Xi) 

a‘^g{ax/X) 


(26) 


Remark 2. The upper bound of the absolute constant used in theorem 8 is uniform over the class Q. 
In specific cases this bound can be considerably sharpened. For example, it is obvious that g{x) = \x\ G Q. 
For such a function g inequality (26) takes the form of the classical Berry-Esseen inequality for Poisson 
random sums, the best current upper bound for the absolute constant in which is given in m- 


Aa 




0.3031 


a^y/X 


(27) 


5 Convergence rate estimates for mixed Poisson random sums 

5.1 General results 

In this section we extend the results of the preceding section to the case where the random number of 
summands has the mixed Poisson distribution. For convenience, in this case we introduce an «infinitely 
large» parameter n G N and consider random variables N* such that for each n G N 

CXD 

P(X* = k) = J e-^^dP{An < A), A: G N U {0}, (28) 

0 

for some positive random variable A„. For simplicity n may be assumed to be the scale parameter of the 
distribution of A„ so that A„ = nA where A is some positive «standard» random variable in the sense, 
say, that EA = 1 (if the latter exists). 
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Assume that for each n G N the random variable N* is independent of the sequence Xi, A 2 ,.... As 
above, let Sm* = Xi + ... + Xv* and if N* = 0, then Sm* = 0. 

From (28) it is easily seen that, if EAjj < 00 , then EX* = EAjj so that DSn = u^EA^. 

Let Nx be the random variable with the Poisson distribution with parameter A independent of 
Xi, X 2 ,... For any x G M we have 

OD 

P{Sn* < xa^/EAn) = ^ P{N* = k)P{Sk < xa^/EAn) = 

k=0 


= J2P{Sk< xax/^) / P(Xa = k)dP{An < A) = 

fc=0 { 

00 00 _ _ 

= I P{Sn, < xa^/^)dP{An < A) = y < x^^yP{An < A) = 


OD 

= y d>( X 

0 


EA^ 

T 


dP{A.n < A) + 


00 

li 


\aVX 


< X 


EA. 


dP{An < X)— ( X 


EA. 


”Vp(An<A). (29) 


A J 


From (29) it follows that 

A* = sup 

a 

00 

^ [ sup 


P{Sn* < xax/EAn) - y ^>(^-^^dP(A„ < AEA„) 
0 

00 


< 


VcrVA 


< X — ^(x) 


dP(A„ < A) ^ / AAdP(A„ < A). 


(30) 


0 0 

Now, if to estimate the integrand A^ in (30) we use theorem 7 and recall the notation F{x) = P(X < 
x), then by the Fubini theorem we arrive at the representation 


A* < 


1.8546 




'EXfmm 


I Vv/Aj 


1-8546 

dP(A„ < A) = — 


00 00 




x^ min 




0 —CO 


1.8546 


00 00 


(7^ 


in |l, J^|dP(A„ < A) 


mm 


ax/X 
dFi{x). 


—00 0 


For X G M introduce the function 


G„ 


xl 1 




a x/K 


I(A„ ^ ^ 
cr^ 


dP(A„ < A) = 

(31) 

(32) 


The expectation in (32) exists since the random variable under the expectation sign is bounded by 1. Of 
course, the particular form of Gn{x) depends on the particular form of the distribution of A„. From (30), 
(31) and (32) we obtain the following statement. 

Theorem 9. If EA„ < 00 , then 


A* < 


1.8546 


cr^ 


EX2g„(Xi) = 


1.8546. 






1.8546 




EXfl{\X,\ ^ ax/K) + E^^l(|Xi| < ax/K) 

'X xj Xfi 


where the random variables Xi and A„ are assumed independent. 

In the subsequent sections we will consider special cases where A„ has the exponential, gamma and 
inverse gamma distributions. 
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5.2 Estimates of the rate of convergence of the distributions of geometric random 
sums to the Laplace law 

In this section we consider sums of a random number of independent random variables in which the 
number of summands N* has the geometric distribution with parameter p = n G N: 


P(K = k) 


1 / n 

n + lVre + 1/ ’ 


fc G NU {0}. 


(33) 


As usual, we assume that for each n G N the random variables N*,Xi,X 2 , ■ ■ ■ are independent. We again 
use the notation Sjsf* = Xi + ... + X]\f*. If N* = 0, then we set Sn* =0. It is easy to see that EN* = n, 
DS'at* = na^. Note that for any A: G N U {0} 


P(iV* = k) = - [ P(A^a = A:) exp I - -}dX, 
n J In) 


where Nx is the random variable with the Poisson distribution with parameter A. This means that for N* 
representation (28) holds with A„ being an exponentially distributed random variable with parameter 
In what follows we will use traditional notation 

OO 2 : 

and TH . r(„.0) = 7(„.oo) 

^ 0 

for upper incomplete gamma-function, lower incomplete gamma-function and gamma-function itself, 
respectively, where a > 0, z > 0. 

In the case under consideration 


OO 

0 0 


where C{x) is the Laplace distribution function corresponding to the density 

£(x) = xGM 

a/2 

(see, e. g., lemma 12.7.1 in [5])- 

At the same time, the function Gn{x) (see (32)) has the form 


r 1 

1 kl [ 

l ncr^ J 

' na J 


jo 


j-A/n 

TT 


X 


dX = 7 1, —^ H--, 


x| _/l x^ 


na^) oJn \2 na 


y 2 ■ 


So, from theorem 9 we obtain the following result. 

Corollary 8. Let N* have the geometrie distribution (33). Then 


Ira/r. a./ M 1-8546 f, 

sup |P(<SAr* < xax/n) — L{x)\ ^- 2 —\ ‘ 

X 






na^ 


+ 


1 


awn L 


|vi7(i.4 
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5.3 Estimates of the rate of convergence of the distributions of negative binomial 
random sums to the variance-gamma law 

The case more general than that considered in the preceding section is the case of negative binomial 
random sums. 

Let r > 0 be an arbitrary number. Assume that representation (28) holds with A„ being a gamma- 
distributed random variable with the density 

y-l -\/n 

pW = rrf ) ^> 0 - 

n^i (r) 

Then the random variable N* has the negative binomial distribution with parameters r and 

OO 

P(N: = k) = | A:ENU{0}. (34) 

^ ^ n'T(r) y k\ T{r)kl \l + nJ \l + nJ ^ ^ 


Let 


Vr{x) = 


4> 


r(r)y VVa 


y ^dX, X G 


be the symmetric variance-gamma distribution with shape parameter r (see, e. g., [28|i. 

In the case under consideration EiV* = EA„ = nr so that DS'tv* = nra'^ and for any x G 


OO _ _ OO 

|4,(.y!^)dP(A„ < A) = ^ /l‘(^\/f. 




r( 

0 

Here the function Gn{x) (see (32)) has the form 
Gn{x) == 


OO 


n^r( 




e-^/'^dX + 


cjre''r(r) 


^r-3/2g-A/n^^ = 


2 ) 2 /(t 2 


1 


r(r) 


X 

7(A —2 
recT^ 




o-Vn 


r r- - 


1 x^ 


2 ’ ncr^ 


So, from theorem 9 we obtain the following result. 

Corollary 9. Let N* have the negative binomial distribution (34). Then 

snp\PIS^,<xa^)-VAx)\ < ^|E[xh(r i)] + ^E[|X.|Y(r - 1 i)] }. 

5.4 Estimates of the rate of convergence of the distributions of Poisson-inverse 
gamma random sums to the Student distribution 

Let r > 1 be an arbitrary number. Assume that representation (28) holds with A„ being an inverse- 
gamma-distributed random variable with parameters | and ^ having the density 


r/2;,-r/2-l 

pW = exp 


2 C2 r(§) 


{-S}' 


A > 0. 
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Then the random variable N* has the so-called Poisson-inverse gamma distribution: 


P(K = k) = 


n 


rl 2 




2 ’-/2r(i 


\ ^ 

-^A"^/^"^exp| - ^jdA, A:ENU{0}, 


( 35 ) 


0 


which is a special case of the so-called Sichel distribution see, e. g., 

n 


. In this case 


EAn = 


so that 


r - 2 


D5* = ^ 
” r- 2 ' 


Nevertheless, we will normalize random sums not by their mean square deviations, but by slightly different 
and asymptotically equivalent quantities a^njr. 

As is known, if A„ has the inverse gamma distribution with parameters | and then A“^ has the 
gamma distribution with the same parameters. Therefore, we have 


n 


r/2 


r(i 


5 _ CXD - - 


2 ^/ 2 r(§) 


d> 


^^‘^dX = Tr{x), 


X E 


where Tr{x) is the Student distribution function with parameter r (r «degrees of freedom») corresponding 
to the density 

r(^) ^ x‘^\-{r+l)/2 


tr{x) = 


1 + 


0 FFr(§) 

see, e. g., d]. 

In this case the function Gn{x) (see (32)) has the form 


, X E IR, 


n 


rjl 


2r/2r(p 


G„(x) = p(a.;‘ > ^) + Me\/aA|i(a;> < A) 


oo a^jx^ 

) J 2-/2ar(§) J 

0-20 


1 

r /r na‘^\ 


'r -\-l na^y 

r(i) 

r V 2 ’2x2J 

^ a \l 2 ^\ 

2 ’ 2 x 2/1 


where 7 (-, •) and r(-, •) are the lower and upper incomplete gamma-functions, respectively. So, from 
theorem 9 we obtain the following result. 

Corollary 10. Let N* have the Poisson-inverse gamma distribution (35). Then 


A* = sup 


P( S'tv* < xa^j-] - Trix) 




1.8546 


2 /r na^\ 


(T2r(§)\^L^^^l2’ 2X2 j 


-E 


1 n 


ay 2 




r -|- 1 
2 
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